CLNS 97/1461 
CLEO 97-2 

Study of the B^ Semileptonic Decay Spectrum at the T(4S) 

Resonance 

(February 7, 2008) 

Abstract 

a^ 

We have made a first measurement of the lepton momentum spectrum in a 

f^i . sample of events enriched in neutral -B's through a partial reconstruction of 

3-^ , i?" — >• D*~£~^i'. This spectrum, measured with 2.38 fb^^ of data collected at 

the T(4S) resonance by the CLEO II detector, is compared directly to the 

inclusive lepton spectrum from all T(4S) events in the same data set. These 

two spectra are consistent with having the same shape above 1.5 GeV/c. 

'T^ . From the two spectra and two other CLEO measurements, we obtain the 

[^^ , B^ and B~^ semileptonic branching fractions, 69 and 5+, their ratio, and the 

^ I production ratio /^ //oo of B^ and B^ pairs at the T(4S). We report b^/bo = 

^ ; 0.950tolll±0mi, bo = (10.78±0.60±0.69)%, and 6+ = (10.25±0.57±0.65)%. 

^^ ■ b+/bQ is equivalent to the ratio of charged to neutral B lifetimes, tj^/tq. 
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It is often argued that the individual semileptonic branching fractions, 6+ and 60, of B~^ and B^ 
should be similar, that a heavy quark meson should decay weakly without regard to the flavor of the 
light spectator quark. Strictly speaking, however, weak hadronic decays are not isospin invariant; 
in the case of charm, the asymmetry of hadronic widths is known to be large. For 6-hadrons, the 
resulting deviation from unity of the lifetime ratio t{B^)/t{B^) is predicted to be of order 10% 
[Q]. There has also been some speculation that differences will be manifested in the semileptonic 
decay spectra jQ]. 

The method of partial reconstruction of B^ —> D*^i^v, where the decay D*~— >Z)'^7r~ is identi- 
fied using only the 7r~, has been used previously by CLEO II, to measure both the mixing parameter 
Xd @Q and semileptonic branching fraction 69 @i^ of B^ ■ We present here an analysis to measure 
the leptonic spectrum of an event sample highly enriched in B^ mesons through the application 
of partial reconstruction as a tag. This is, with minor corrections, a measurement of the inclusive 
leptonic spectrum from semileptonic B^ decays. It may be compared with the inclusive spectrum 
in T(4S) events, which contains a roughly equal mix of B^ and B^ . By comparing the rates in 
untagged and tagged events, it is possible to extract not only a value for 60, but also for 6+ and 
(5 = 6+/60. Assuming that the exclusive and inclusive semileptonic partial widths are equal for B^ 
and i?"*", (3 is then equal to the ratio of charged to neutral B lifetimes Tj^/tq. 

The data, collected with the CLEO II detector p| at the Cornell Electron Storage Ring (CESR), 
consist of integrated luminosities 2.38 fb~^ on the T(4S') resonance and 1.13 fb~^ at a center-of- 
mass energy which is lower by 60 MeV (continuum). All events considered here must have at least 
5 well-fitted charged tracks, a measured energy at least 0.15 times the CM energy and an event 
vertex consistent with the known interaction point. 

B^ decays are identified in the mode B^^D*~'i^v [D*" —>■ D^tt^) (inclusion of the charge 
conjugate decay is implied throughout this paper), where only the £+ and vr" are detected. A large 
sample of B^ is obtained by exploiting the extremely low energy of the D*~ decay. The method is 
described fully in Ref. Q. In brief, the momentum of the it~ is scaled to obtain an approximate 
four-momentum {E£)*,pu*) for the D*^: 

Ed' ^ ^Md* = Ed* 



PD* =^ Ptt X \/E%, - M^, = PD* 

where E-,^ is the vr" energy, E^ = 145 MeV is the energy of the vr" in the D*~ rest frame, and 
Mjj* is the mass of the D*~ . Using the approximation pB — allows us to calculate a squared 
missing mass, 

^' = (^beam " ^^* " ^^f " (P^* + P^)' 

which approximates the squared mass of the neutrino. Nearly all £~^tt~ candidates stemming from 
decays B -^ D*^ {■k)1'^ v (D*^ -^ D^n^), where the pion is the daughter of D*~ , fall in the signal 
region, M^ > —2.0 GeV^. Such candidates constitute our tag while all other £^tt^ combinations 
constitute the background. 

The number of tags is counted in the entire data set and in subsets which contain additional 
leptons. To avoid confusion, the lepton used in the tag candidate will be referred to hereafter as the 
"tag lepton" , while the additional lepton from the remainder of the event will be called the "second 
lepton". We allow more than one "second lepton" to exist in an event. By sorting event subsets in 



100 MeV/c increments of the momentum of the second electron or muon, we measure their spectral 
distributions, which may then be compared to the corresponding inclusive lepton spectra 0. 

Each tag lepton track must have a momentum in the range 1.4-2.4 GeV/c and be consistent with 
originating at the interaction point. Electrons are identified primarily by the ratio of calorimetric 
energy to momentum and specific ionization {dE/dx) in the drift chamber. Muons must penetrate 
at least five nuclear absorption lengths. 

The pion in a tag candidate must be consistent with originating at the interaction point, have 
charge opposite to that of the lepton, and have momentum less than 0.190 GeV/c. If dE/dx 
information is available, it must be consistent with the pion hypothesis. 

Second electrons (muons) are required to have momentum between 0.6 GeV/c (1.5 GeV/c) and 
2.5 GeV/c. Since B^ mixes with B^ , we accept second leptons of either the same or opposite sign 
relative to the tag lepton. The tracking criteria are identical to those for the tag lepton except 
that the track polar angles are more restricted, jcos^^l < 0.61 for muons and jcos^^j < 0.707 for 
electrons, where 0^ is the track polar angle. By matching the tracking and identification criteria to 
those used in the determination of the inclusive spectra 0, we minimize systematic uncertainties in 
our comparison of the two. To eliminate multiply reconstructed tracks, a second lepton candidate 
is rejected if the angle between its momentum and that of the tag lepton is small, cosOm > 0.99. 

To obtain the net number of tags in each event subset, we select candidates in the signal 
region (M^ > — 2.0 GeV'^) and subtract backgrounds. The nonresonant contribution is estimated 
using the continuum dataset, corrected for luminosity and energy differences. Where a second 
lepton is selected, we account for cases where this second lepton is a misidentified hadron (fake) by 
repeating the above process selecting, in place of second leptons, tracks not identified as leptons but 
passing all other second lepton requirements. The result is scaled by the fake rate per track at the 
appropriate momentum 0. After correction for fakes, the remaining background, from candidates 
which are random combinations of leptons and pions in BB events, is estimated via Monte Carlo. 
Its normalization is obtained from the sideband region, —20 < M^ < —4 GeV^. The set of events 
used for this estimate is approximately five times the size of the data set. The bins with second 
lepton momentum in the range 2.1-2.5 GeV/c are combined for the Monte Carlo sample to improve 
the statistical significance of the background estimate. The statistical error of the Monte Carlo 
sample is included with that of the data and contributes approximately 10% to its magnitude. Tags 
that are counted more than once within the signal region because the soft pion curls and is multiply 
reconstructed are found in Monte Carlo to occur for about 4% of the signal. This overcounting 
cancels in the ratio of tagging rates which is taken for the final result. The resulting overestimate 
of the statistics is accounted for by increasing the statistical errors by 2%. 

The evaluation of background in the signal region depends on an accurate shape for the simu- 
lated background distribution in M^. To quantify the uncertainty on this, we repeat the analysis 
using "wrong-sign" candidates, where the tag lepton and pion carry the same charge and no net 
signal is expected. We obtain a net "signal" of —277 =b 211 events from the entire event set. This 
is consistent with zero, and we take as the systematic uncertainty the absolute value of the mean 
plus one standard deviation, which corresponds to 3.0% of raw candidates in the wrong sign signal 
region. We then take the uncertainty in the right sign background of each event subset to be 
the same fraction. The resulting systematic errors for tags with a second lepton are found to be 
1.4%-2.1% of the net signals. We assign an overall error of 2.1% to the second lepton spectrum. 

The net number of tags in the entire data set, Nt, is 22, 880 ± 283 ± 491. Figure | shows M^ 
distributions with estimated backgrounds for the full data set and for subsets corresponding to 



several lepton momentum bins. The net number of tags with second leptons above 1.5 GeV/c is 
1268 =t 54 (statistical error only). 

The net number of tags in each bin of second lepton momentum is now corrected for the 
lepton identification efficiency. We use the combined efficiency for geometric acceptance, track 
reconstruction, and identification as a function of momentum which is found in Ref. 0. 

We now estimate and subtract the contribution from ^ decay. The spectrum of electrons or 
muons from ■(/; decay which accompany selected tags is found using Monte Carlo BB events. It is 
then scaled to the number of tp^s in data. The scaling factor is determined through reconstruction 
of decays '0 — > fJ''^lJ'~ in the full event set and in Monte Carlo. 

The remaining second leptons can be considered as either primary, produced directly in semilep- 
tonic B decay, or secondary, from products of B decay. To assess the contribution from secondary 
decay, the tagged electron and muon spectra are now fitted simultaneously to sums of primary and 
secondary shapes, as was done for the inclusive spectra 0. The primary shape is derived using 
the model of Altarelli et al. S. Secondary leptons originate from decays of charm (Z?, Dg, Ac) and 
T, as well as (for electrons) n^ Dalitz decays and conversions of photons in the beam pipe. Above 
0.6 GeV/c, the spectra from all of these are similar in shape. They are therefore approximated 
using the spectrum from D's. The inclusive spectrum of Ref. Q, which was explicitly corrected 
for all but the D contribution, was fitted using the same shape. The tagged spectra, normalized 
to Nt, are shown with fits in Figure ^. Their overall appearance differs from that of the inclusive 
spectra because the secondary rate is higher relative to the primary rate; as one primary lepton is 
in the tag, there is at most one other primary lepton, whereas there can be two secondaries from 
charm. 

The secondary contribution is now subtracted from both tagged and inclusive spectra. We then 
apply two corrections to the tagged spectrum. The first, arising from the efficiency of requirements 
on cosOu and on track multiplicity, has a net value of +1.1 it 1.1% [H. The second arises because 
the tagging efficiency has a slight dependence on the decay charge multiplicity of the other B in 
an event, such that if the other decay is semileptonic, its lower average multiplicity results in a 
higher tag efficiency. The correction was calculated for each momentum bin of the second lepton 
using a Monte Carlo simulation and found to be independent of second lepton momentum, with 
an average value of (—2.9 it 1.3)%. The net correction to the tagged rate from the two sources is 
(-1.8 ±1.7)%. 

The tagged spectra of electrons and muons are then divided by the corresponding inclusive 
spectra, bin-by-bin. If the primary spectra have the same shape, the resulting distribution should 
be flat, with a value equal to R = btag/binch where btag and bind are the semileptonic branching 
fractions measured in the tagged and total events samples, respectively. The distributions appear 
indeed to be flat. They are thus fitted to a constant function, where we now restrict the momentum 
range to 1.5-2.4 GeV/c in order to limit uncertainties associated with the secondary spectrum 
(Figure |). The fits yield R = 1.013 ± 0.061 (1.036 ± 0.071) with x^/DOF = 7.5/8 (3.1/8) for 
electrons (muons). The weighted average is 

R = 1.023 ± 0.046. 

We also obtain R by taking the ratio of the primary contributions to the spectra from the fits 
described above, and find excellent agreement. 

Because the procedures for obtaining tagged and inclusive spectra have been made nearly 
identical, many of the associated systematic errors are highly correlated and cancel approximately in 
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the ratio R. Among these are errors from lepton tracking and identification efficiencies, subtraction 
of the if) contribution, fake rates, and spectral shapes for fitting to find the secondary contribution. 
To estimate the effects of variations in the ip contribution and fake efficiencies, we repeat our 
measurement, adjusting each contribution in turn to the hmits of the estimated uncertainty. To 
assess the uncertainty due to shapes used in fitting to obtain the secondary contribution, fitting was 
repeated with varying primary and secondary shapes. For the primary shape we include the model 
of Isgur et al. 0] where the fractional contribution of i? — > D**lv to semileptonic decays is fixed 
for each fit in the range (23 it 1)% [^ (ISGW**). For the secondary spectrum we use the method 
and shapes used in Ref. H . These correlated sources and effects are summarized in Table |. 

Uncertainties which apply to the tagged but not to the inclusive measurement include effects due 
to event selection, angular cuts, and tagging efficiency differences discussed above, which contribute 
a total of 1.7%. The uncertainty from the background estimate for the tag cancels approximately 
in dividing by Nt, and we assign an overall error of 1% to the magnitude of the tagged spectrum 
from this source. In addition, we approximate the spectra of all secondary leptons in the tagged 
spectra by those from D decay only. To estimate the uncertainty from this and from the fact that 
there may be small differences in shape between tagged and inclusive samples due to tag selection 
itself, we repeat the analysis after subtracting an estimated contribution of secondaries from sources 
other than D decay. This causes a shift of 1.7%, which is taken to be the error. The net systematic 
error on the value of i? is 3.5%. 

R may be represented by an equation in terms of ftoi ^+) 13 = b+fbo, and several other parame- 
ters, which can be solved for bo, 6+, and f3. The average semileptonic branching fractions of B's in 
the tagged and inclusive event samples, btag and bind, may be expressed as 

bind = foobo + f+-b+ = foobo{l + a) (1) 

btag = fo,tagbo + f+MgK = &o(l + f+,ta9^P), (2) 

where /oo and f^ are the fractions of B^ and i?^ events from the T(4S), fo.tag and f+,tag ^^^ ^^^ 

corresponding fractions among tags, a = /-| 6+//00&O) and 5(3 = j3 — 1. 

From (|2|) and the definition of /3 we obtain 

J '^to.g ^ "ind /Q\ 

° " 1 + f+,tag6P ~ 1 + U,tag6P ^^ 

Rb,nd{l + Sp) 

To obtain 5/3, we take the ratio of (^) and (|l|) 

~ bind /oo^o(l + a) 

_ 1 + f+,ta9^P , . 

^^°°- R{l + a) ■ ^^^ 



From the definitions of a and 5/3, 



Equating (^ and (0), 



_ 1 + 5/3 



(1 + 5R){1 + a)(l + 5P) = (1 + U,tagSP)ia + l + 5p) (8) 

f+^tag, SP, and SR = R — 1 are expected to be "small", < 0.1, so we discard the highest order term, 
f+^tagSP'^, to arrive at a linear equation in SP which is solved to obtain 

For a 

a = B{T{AS) -^ 

M 

B{T{4S) -^ B+B- 

get a = 1.15 ±0.13. 

The fraction f+,tag i^ nonzero because some charged i?'s are present among our tags from 

decays of the type B -^ D*^'Ki^v (where the D*^tt may or may not form a resonance). The 

fraction /** of such decays in the tag sample may be determined by adding the reconstruction 

of the D^ decay, so that a measured D*~ momentum is used in the calculation of missing mass 

squared. In this case the resolution is sufficient to achieve some degree of separation between the 

decays without and with an additional pion, so that the relative rates may be determined. In 

terms of our analysis the only directly applicable measurement has been performed at CLEO I, 

with the same event and tag lepton momentum requirements which found /** = 0.14 it 0.08, the 



"" {l+a)U+M9-5R)-<^' 


K-^J 


we average the CLEO 


II measurements 


B+B-,B+ -^ D*^^+v)/B{T{AS) -^ B^B^,B^ -^ D* 


-i+u) = 1.14 ±0.14 ±0.13 


and a 


= 


^B+ -^ ipK^*^+)/B{T{4S) -^ BW^,B^ -^ V^WO) = 


1.15 ±0.17 ±0.06 [11|, to 



value of which can be derived from measurements cited in |12|. Because /** is a ratio, it depends 
mainly on kinematic criteria and should not differ significantly in CLEO II. We use a measurement 
from ALEPH, B{B -^ D*-Tr+£+uX) = (1.25 ± 0.16 ± 0.12)% H, and Monte Carlo simulations to 
obtain an independent value, /** = 0.10 ± 0.02, which is in agreement with CLEO I but may have 
a large systematic error due to insufficiency of our current knowledge about the different modes 
which contribute to S — > D*~'Tr^i~^L'X. We use the CLEO I value because it is more directly 
obtained, and our evaluations of bo, 5+, and /oo are insensitive to it. The relative contributions to 
/** from charged and neutral -B's is determined using the approximate conservation of isospin in 
the fragmentation of the hadronic part of the final state; if 6o and 6+ are equal and the charged and 
neutral S's decay to D*Xi^u at the same rates, the rate to states where X is a charged pion is 
expected to be double that of the rate where it is a neutral pion. In other words, B^ — > D*^'k'^I!.^u 
contributes to /** at twice the rate of B^ -^ D*~TT^i~^i', weighted by the respective production 
rates and semileptonic branching fractions 0, 

2a 
l+^a' 



f+,taa = r-r-iTZ.^ (10) 



from which we get f+,tag ~ (^-^ ^ 5.5)%. 

Using R = 1.023 ± 0.046 ± 0.036, /** = 0.14 ± 0.08, and a = 1.15 ± 0.13, we obtain 

a -n qc;n+0.117+0.091 

/+_//oo = 1.21 ±0.12 ±0.17, 

where in each case the first error is statistical and the second systematic. The latter number may be 
expressed /oo = 0.452 ± 0.028 ± 0.038. The values of ho and 5+ may be obtained with the addition 
of hind- The most recent published measurement of an inclusive semileptonic branching fraction at 
the T(4S) was made using a lepton tagging method which gives a value that is not strictly equal to 

8 



bind [141, as the tagged event sample may not have the same composition as generic T(4S) events. 
However, given known hmits on the input parameters, it is expected to differ from bind by less than 
one part in 10^, so we use the result, bind ~ (10.49 ± 0.46)%, to get 

6o = (10.78 ± 0.60 ± 0.69)%, 
b+ = (10.25 ±0.57 ±0.65)%. 

As the determinations of bQ, b+, j3, and /-| //oo involve correlated inputs, the errors are obtained 

numerically. All of the statistical uncertainties are defined by the uncertainty in R. The dominant 
contributors to the systematic uncertainties are R for /5, both R and bind for 6o and 6+, and a for 

/oo- 

Our results are consistent with our previous ones 1^, /? = 0.93 ± 0.18 ± 0.12, 6o = (10.9 ± 0.7 ± 
1.1)%, and 6+ = (10.1 ±1.8± 1.5)%. Our new result for 6o is somewhat correlated with the previous 
one, as the data sets overlap partially and parts of the method are similar. The results for 6+ are 

uncorrelated. /_) //oo may also be obtained using the CLEO measurement of a with a lifetime 

ratio based on measurements at LEP and CDF, as in [^. The value reported here was obtained 
using CLEO data exclusively. 

To summarize, the electronic and muonic spectra have been measured in a sample of events 
enriched in neutral i3's using partial reconstruction tags. By comparing these spectra to the 
inclusive spectra in T(4S) events, we probe differences between B^ and B^ in rates and spectral 
shapes of primary semileptonic decay. We observe no evidence for differences in the spectra between 

1.5 and 2.4 GeV/c. Assuming the shapes are the same, we obtain 6o) b+i Pi and /-| //oo- Our new 

value of /-I //oo is the first that is independent of B meson lifetime measurements. 
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TABLES 

TABLE L Systematic errors on R. Those listed in the upper section apply both to tagged 
and inclusive spectra, while those in the lower section apply only to the tagged spectra. For each 
listed source, the ratio of its magnitude to the net tagged spectrum above 1.5 GeV/c is shown in 
parentheses. 

Source asys/value{%) 6{R)/R{%) 

or procedure 

i/j contribution (0.011) 30 (L3 

Fake electrons (0.002) 50 0.1 

Fake muons (0.011) 25 0.2 

Secondary spectrum (.044) as in Ref. Q 2.3 

Primary spectrum (0.95) ISGW** 0.3 

(6 -^ u) (0.019) 21 - 

Event/tag selection 1.7 1.7 

Background estimate (0.59) 1.0 1.0 

Secondary spectrum (0.044) vary sources (see text) 1.7 

Net systematic error 3.5 
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FIG. 1. Tag candidate M^ distributions, data (symbols) and simulated backgrounds (solid 
histograms). Continuum has been subtracted in all cases, and in (c)-(f), fakes have been subtracted. 
(a)-(b): full event set (errors are smaller than symbols), (a) opposite-sign candidates, (b) same-sign. 
(c)-(f): events containing a second lepton, opposite-sign candidates, selected bins as indicated. 



12 



3521296-002 



0.04 



^^ 0.02 



> 

(D 

- 

^ 0.04 

Q. 

■a 



0.02 



— 






w — 




— 1 \ \ \ \ 1 \ \ \ \ \ 

Electrons 

'V^ 


1 1 1 1 1 1 1 1 1 1 l'^ ^1 " 1 1 1 


— 








1 


1 1 1 1 1 1 1 1 1 

Muons 


1 1 1 1 1 1 1 1 1 1 1 ^l'-' 1 1 1 



12 3 

Momentum (GeV / c) 

FIG. 2. Measured spectrum of second leptons, corrected for efficiency and with ip contribution 
subtracted, shown with fit to primary (dashed) and secondary (dotted) spectra. 
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FIG. 3. Bin-by-bin ratio of tagged and inclusive spectra, with fits to a constant (dashed line) 
for electrons (top) and muons (bottom). 
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